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I. INTRODUCTION 



It is well known that particles immersed in a thermal medium at temperature T have 
their properties modified. The knowledge of the reasons for these modifications and their 
consequences is of great relevance in the context of finite temperature field theory (FTFT). 
However, quantum field theory at finite temperature has some subtleties such as the break- 
down of the perturbative expansion which is manifested in two cases. The first one is the 
appearance of infrared divergences in massless field theories or tachyonic masses in theories 
with spontaneous symmetry breaking [Q-Q]. The second breakdown surges in calculations 
using higher order diagrams where large T compensates for the powers of the coupling con- 
stant. The practical way to circumvent these problems is a resummation of certain classes 
of graphs up to a given order in the perturbative expansion. This resummation (of an 
infinite set of diagrams) has been studied intensively in the last three decades 
Nevertheless, the use of these kind of nonperturbative methods does not altogether solve 
the problem of the ultraviolet (UV) divergences (which will become temperature dependent 
through the gap equation for the mass) that are still present. In the related literature, it is 
frequently claimed that the poles are suppressed by the renormalization procedure, but it is 
not always shown how. It is not necessary to say that, if one wants believable results, the 
problems raised by the divergences must be satisfactorily resolved. For example, consider 
the contribution to the mass squared at one loop in the neutral model, 

where mo is the vacuum mass, Un are the Matsubara frequencies, defined as Un = 2mTT for 
bosons, nB^x) = l/{e^^ — 1) is the usual Bose-Einstein distribution and ujk{m), the rela- 
tivistic energy, is given by \/k^~^^m?. The first term in the integral is the finite temperature 
contribution, which is free from UV divergences. The last term in the integral is the zero 
temperature part (at this stage) and is divergent. By the use of some non-perturbative 
method like the Hartree approximation, one gets M^(T) = rriQ + Ili{M), with M being the 



2 



dressed mass. If the mass running in the loop is the vacuum mass (i.e., as in conventional 
perturbation theory), then the theory can also be renormalized at finite temperature Q but 
if some resummation approximation is naively used as in the example above, the theory is 
non-renormalizable at finite temperature since one can not absorb in the original Lagrangian 
the infinity term, A ^^^'p | (obtained using dimensional regularization [|TU|)- The nonobser- 
vance of this "incompatibility" is one of the reasons for the failure in finding appropriate 
renormalization conditions for a theory in study. The difficulty in obtaining renormalization 
prescriptions in a theory treated by a self-consistent approximation method at finite tem- 



perature is also manifested in other theories such as the linear sigma model [|lT|-[l^. Being 



aware of these problems, the authors of proposed to treat the linear sigma model as 



an effective model, with the undesired result that the theory and its predictions are sensitive 
to the value of the finite momentum cutoff A. In fact, it was known more than 20 years 
ago that these approximations requires temperature-dependent counterterms [0 . Then one 
concludes that in order to achieve renormalization in self-consistent approximation methods, 
the counterterms need also to receive the benefits of resummation JIB|- p!B| , pT| . Of course it 



is not necessary to resum the counterterms if one employs ordinary perturbation theory at 
finite temperature to leading order in the perturbative expansion [0,^ since finite temper- 
ature does not introduce any new UV divergence to the theory [0. In a recent paper ||21 



we presented a Modified Self-Consistent Resummation (MSCR) which resums higher-order 
terms in a nonperturbative way and cures the problem of breakdown of the perturbative 
expansion. Since this method is developed in steps, it has the essential features that allow 
the absorption of UV divergences, the avoidance of overcounting of diagrams and the identi- 
fication of regions of validity |21[]. The main difference between the methods used in [p!^ -[T8[] 



and MSCR are that in W& and [IS] a mass parameter was introduced in the beginning to be 



^Considering that the theory under investigation is renormalizable at zero temperature. In the 
case of the Xcj)'^ model one has to consider also the "triviality" behavior of the theory p] . 
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determined later by some criteria and in [|T^] a known mass parameter was added (having in 
both cases the subtracted term treated as a new two-point interaction), whereas MSCR keep 
the same fundamental theory in the recalculation of the self-energy |^. It is important to 
point out that, although one has, in principle, the freedom of adding and subtracting mass 
parameters to a LagrangianQ one must be careful. For instance, in the case of (where 
a resummation was applied in the linear sigma model at one-loop order in the perturba- 
tive expansion), this operation provided the pion and the sigma dressed masses with the 
same corrections. This is incorrect, since the self-energy for these fields are different, at 
least at one-loop order. It is worth to remark that even if they had added different (and 
correct) contributions to the sigma and pion fields, neither the 0(4) linear a model would 
be renormalizable nor Goldstone's theorem would be satisfied in the region of intermediate 
temperatures (around Tc, the critical temperature) This is because in that temperature 
region, quantum fluctuation may need a more detailed description and certainly, more than 
one-loop in the perturbative expansion |2^. In this paper we investigate the neutral scalar 
model at finite temperature. We apply MSCR, and use imaginary time formalism to 
calculate the self-energy up to two-loop order in the perturbative expansion. Some previous 
works describing two- loop calculations may be found in P^ , p3[| and more recently in [2^. We 
follow the conventions of Ref. The reader who is unfamiliar with FTFT may consult 
and the more recent publications |]25|j26| . 

This paper is organized as follows. In Section II we calculate the self-energy of the 
model up to two-loop order and isolate the divergent and finite temperature contributions 
in order to proceed with the resummation and renormalization. In Section III we briefly 
note MSCR which will execute this resummation. We conclude in Section IV. 



^In the classical level this procedure changes nothing. However, to finite order in perturbation 
theory it has consequences 
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II. THE LABORATORY MODEL 



Consider the following Lagrangian density describing the self-interacting A0^-theory: 

C = C,{<P) + C'^\<P) + C^\(t^,X), (2) 

where 

= ^d.ct^d^ct^ - — </)^ (3) 
the interaction Lagrangian is given by 

= (4) 

and the counterterm Lagrangian £'^*(</), A) is expressed as 

A) = \A{X)d,cj>d>^ct^ - \B{X)m^^ - ^C(A)0^ (5) 



As we shall see next, up to two- loop level the coefficient of the counterterms are [0,|T7 



~ (47r)2 2e ^ (47r)n2e2 4g 

_ A' 3 A^ / 9 3 
~ (47r)2 2e ^ (47r)n4e2 2e 



where we use the modified minimal subtraction {MS) scheme. Through out this paper, 
we employ dimensional regularization, but we omit, for notational simplicity the factor /i^*^ 
which multiplies A. In the Xcj)^ model, only the self-energy has to be resumed P,|17|- So, 



the coupling constant A can be treated as the expansion parameter without the necessity of 
resummation. This allow us to make the calculations in terms of the renormalized A. 

Let us now compute the one and two-loop contribution to the self-energy of 0. The only 
contribution to the self-energy at one-loop comes from the diagram in Fig. 1 and it is given 
by Eq.(|l]). One obtains the self-energy in the usual fashion 0, 
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A ml 
2(47r)2 



-- - 1 + ln 

e 



ni(mo) = n?(mo) + nf(mo) = 
A dpp'^UBiujp) 



(7) 



rrin 



+ 



2 Jo 27r2 a;p(mo) 



whith i = — 7 + ln(47r), where 7 is the Euler constant and is the renormahzation scale. 
The divergent contribution of this graph comes only from 115(^0) and will be difined as 



nr(mo) 



X ml 1 



2(47r) 



(8) 



while the finite temperature contribution is given by 



n^(mo) 



A r°° dp p"^ nsiuj) 



2 JO 27r2 c<jp(mo) 



(9) 



The mass and vertex counterterm generated two extra diagrams to one-loop order |T0 
which are drawn in Fig. 2. These diagrams, defined as n^^'*'"'^^ and n^^'*^'"^, will be necessary 
to cancel the infinities formed with the mixing of T = and T 7^ pieces of Fig. 3 and are 
easily calculated to give: 



extra 



X^ml 1 



4(47r)2 g 



(4vr) 



+ 



1 d 



2mo dmo ^ ^ 



(10) 



with 



extra^ 



3A^ ml 
' 4(47r)4 



~2 g 



mmo)=fi''J 



dq 



,3-2e 



3A2 1 



4(47r) 



(11) 



(12) 



(27r)3-2^o;g(mo)' 

At two-loop order there are two contributions 112,1 (mo) and 112,2(^0, K) that are depicted in 
Fig. 3. We note here that only the second two-loop contribution depends on = {k^, k), 
the external four-momentum. The "double scoop" diagram is 



n2,i(mo) = (mo) + i(mo) + (mo) 



(13) 



where 



(mo) 



A^ ml 
4(47r)4 



finite terms, 



(14) 



J d_ 

2mQ drriQ 



rrin 



(/|(mo))+ln ^2 



(15) 



A2 1 



4(47r)2 g 



(16) 



The "setting sun" diagram gives the following contributions: 



ReU2,2{mo, K) = 2(mo, K) + 0^,2(^0, K) + njf (mo, iT) 



(17) 



where we have defined 1122(^05-^); 1^22(^0)-^) and n2'2 (^tiq, -f^) as the real parts of the 
functions Gq{Kq, K)^ G2{Ko,K) and Gi{Kq,K) respectively since, accordingly with the 
notation of |T^, 112,2 ('^o; K) = G0 + G1 + G2 ■ These real retarded parts have been obtained 
after the necessary analytic continuation ko —>■ —ifl + e, with e — >^ 0^ and > 0. Now, with 



Ultimo, K^) = ^^^ 

2,2V 0, ; ^^^^y 



(47r)4 24e 



z + finite terms^ (18) 



Ill2{mo,in,k = 0) 



dp p 



nsiujp) 



(27r)4 Jo ujp{mo) Jo ujq{mo) 



aq q — : r m 



Y_ 



where 



(19) 



X [Q"^ - {Ujp -UJq+ Wpig)^] [fi^ - {Uq + UJp- UJp±qf], 



(20) 



n°|(mo,2fi,A; = 0) 



/|(mo) 
2(47r)2 



1 , ml\ 



+ F(^]^ 



where F(f2^) is 

F(^]2) 

with 



(27r)4 Jo o;fc(mo) Jo ujq{mo) 



dq 



gin 



X4 



-4A; 



(21) 



(22) 



X±(ll ) = [f2 - {uJp + UJq + UJp±q) ] \VL - {Ujq - Up + UJp±q) ]. 



(23) 
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In Eqs. (|I6|) and (plD, Ilg^f (mo) and 112^2(^0) contain the mixed temperature-dependent 
divergent contributions from the T = part from one loop with the T 7^ from the second 
loop in each graph of Fig. 3. When the extra diagrams are taken into account, these spurious 
(not absorbable by counterterms) temperature-dependent divergences fortunately cancel as 
the reader can easily check combining Eqs. (|I6|) and ( |2TD with equation (p!T|). In addition, 
the reader can see that the term In (^^^ in the simple pole has been also canceled |^ . 

Thus, the divergent contribution to the mass at two-loop order after the cancelation of 
the divergences from the extra diagrams with the ones from 112,1 (mo) and 112,2 ('^o, -^)5 is 
given by 

The term j^^^ in Eq.(|l^) is canceled by the usual two-loop wave function renormal- 
ization counterterm whose coefficient is given by A{X) in Eq.(^. The complete divergent 
contribution up to two-loop order to the squared mass is then expressed as 

n'^™(mo) = nf (mo) + nf (mo) = (25) 
A mo 



which justifies -B(A) in Eq.(|^). An important remark is now in order. As we shall see in 
the next section, with MSCR the counterterms that remove divergences in the self-energy 
at zero temperature also remove the divergences at finite temperature. The only difference 
is that the recalculation of the self-energy changes the vacuum mass by a thermal mass. 
However, the cancellation of the temperature-dependent infinities is still guaranteed due to 
the consistence of the MSCR method. 

The temperature-dependent part of the self-energy at one and two-loop order at zero 
three-momentum is given by 

n'^(mo, k = 0) = n^(mo) + nf,i(mo) + Ultimo, ifl, k = 0) - (26) 



2' 



2(47r)2 

where the two terms in the second line of Eq.(^) come from Il2'2{mo,iQ, k = 0). 
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III. DRESSING THE MASS 



We now begin our calculation of the dressed mass up to two-loop order in the perturbative 
expansion. The thermal mass is defined to be the real part of the pole of the corrected 
propagator at zero (k) momentum 

D-i = I)-! + n(mo, zfi, A? = 0) = ^ (27) 
Q'^ = ml + n(mo, ifl, k = 0), 

where T>q is the tree-level propagator and 11 = -1-11^ is the self-energy at finite temperature 
up to a given number of loops at zero momentum. In this work 11° is given by Eg . (p5D and 



is given by Eg . (p6|) . Namely, the finite terms of Il'j', Ilg ^ and 112 2 been dropped. 

However, as we have discussed earlier, the perturbative expansion breaks down in the 
temperature dependent part of the self-energy mainly in two ways. The necessity of a 
consistent resummation is evident both at low and high temperatures. Of course this break- 
down is also manifested in the model as well. The first necessity surges when there 
is a symmetry breaking, and in this situation we would have for the shifted- model: 
rriQ rn? = mg + |Az/^(T), where z^(T) is the thermal expectation value of defined in 



T = as \J with mg < to allow symmetry breaking. In this case, the mass running 

in the loops become tachyonic even below since z^(T) decreases as T increases. The second 
necessity is due to the fact that higher order diagrams are larger than the lower ones at high 
T, even if the strength of the coupling is small. For instance, consider the diagram 112 i (mo). 
In the high temperature limit, it gives A^^. For a diagram with j-bubbles attached, the 
power of ^ become more severe, as A-^ ■ Thus, higher order terms must be resummed 
to get sensible results at high T. 

With MSCR, the resummation is consistently achieved by the recalculation of the self- 
energy. Let us now review our procedure which resums higher loop diagrams in the tree-level 
propagators. The goal is to make renormalization possible when obtaining the pole of the 
effective propagator. The method consist in recalculating the self-energy, in steps, using in 
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each step the mass obtained in the previous one such that = (A„ + 1)M^„^ + n(M„„i), 
where n is the order of the nonperturbative correction and An is the coefficient of the 
appropriate counterterm. With this procedure it is easier to identify and absorb the divergent 
parts of the self-energy (since the masses which multiply the divergences are necessarily the 
same as in counterterms) in order to have finite thermal masses. 
Step 1: 

Start with the effective Lagrangian where the mass is given by: 

Ml = ml (28) 

Step 2: 

Evaluate the one and two-loop self-energy corrections to this mass from Eqs.(p5D and 
(^) and define the ffist order corrected mass as 

Ml = Ml + n(mo, fi, = 0) = {Ai + l)m^ + n(mo, Mq, = 0) = (29) 

m^ + n^'="(mo,Mo,fc = 0), 

where Ai = + - ^) is given by the second line of Eq.(D. 

Step 3: 

Now we take the mass computed in the previous step and improve the results by defining 
a next-order nonperturbative correction. With this we get a new effective Lagrangian where 
the mass is given by 

, " ^ 

Ml = Ml + n(Mi, VL,k = {]) = ml + n^^"(mo, Mq, A? = 0) +n(Mi, = Mi, = 0) = (30) 
{A2 + l)ml + {B2 + l)n^^"(mo, Mo, ^ = 0) + n(Mi, = Mi, = 0) = 

ml + n^""(Mi, n = Mi,k = 0), 

where A2 = Ai. The coefficient of the temperature dependent mass counterterm B2 is fixed 
in a manner to cancel not only the divergence proportional to n(mo, Q,k = 0), but also this 
term together. Generalizing, at each stage of the procedure, for n > 1, in the expressions 
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for M„, the self-energy n(M„_2) have to be cancelled to avoid the overcounting of diagrams. 
This implies that B2 = A2 — I. It is shown that this first recalculation (iteration) corresponds 
to a daisy sum P,[^,pl[] . 
Step 4: 

Proceeding with the iteration, in the limit n — 00 the mass M„ have formally the same 
expression as the mass M„_i which is already renormalized. Thus, in this limit we will have. 



Mi 



m^ + n^^"(M„,fc = 0). 



(31) 



At each intermediate step, in the loops we set Vt^ = M^_i in the computation of M^. This 
ensures the cancellation of the divergences in all stages of the process, since the masses in the 
counterterms will necessarily be the same as in the divergences. In the end, in the resulting 
integral equation of interest, = as it should. A part of the complete self-energy 
diagrams in the superdaisy sum is shown in Fig. 4. As we are concerned only about finite 
temperature effects, n^^"(M„, A; = 0) is given solely by Eq. Strictly speaking, in Eq. 
(^) should be obtained numerically since 11'^^"' (M) cannot be evaluated in a closed form. 
However, we can get a rough expression for the gap equation in the high temperature limit 



12(47r) 



1 - 



AT2 








24 




3M\ 




7tT ) 


^ 24(47r)2 



+ 



327rM \ 12 Atc 



(32) 



In — + 5.0669... 



where the last term in the equation above is the computation of Ilg 2(M) + F{M) which we 
borrow from [|r7|. Thus, the resummed thermal mass up to two- loop is the positive root of 
the equation 



"^° + ^ + 24(4^ ^HtJ 



2A ^ 



A^ 

3847r 



0. 



(33) 



We have shown that |^ to one loop-order the thermal mass of the massless Xcj)'^ in the 
weak coupling limit (A << 1) at high temperature can be dressed algebraically. This simple 
application exemplify that, indeed, the method works very well. The main results are quoted 
here: 
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= Ml + n(Mo 



Mo = 0, 
AT2 



(34) 



24 



M| = M( + n(Mi) = {A2 + i)n^-^"(Mo) + n(Mi 

3Afi 



(35) 



AT2 
24 



AT2 
~2A 



ttT 
A ^ 



247r2 



The mass dressed by the second iteration, M2, which was obtained in our method evaluating 
Fig. 1 with Ml in that loop can equivalently be achieved summing an infinite set of "daisy" 
diagrams with Mq in the loops. In this case all "daisy" diagrams are IR-divergent since 
Mq = 0, but their sum is IR finite [^|17|,^. Continuing the iterations, 



AT2 
~2A 



Ml = Ml + n(M2 
3Mi / 3Mi 



(36) 



AT2 
~2A 



TlT 

A ^ 



1 



247r2 



+ 



ttT 
9 / A 



2 \ 247r2 



and finally, for the ra-th iteration, we obtain 



24 



1 + Ej=l 



247r^ 



The 



■'superdaisy" sum (which corresponds to the limit n — > 00) gives 



M^ 



AT2 
~2A 



A 



2 \ 247r2 



(37) 



1 + 3 (^\ 
/ + 2 V247r2 ) 

As one can see, the nonperturbative nature of the procedure leaves its signature in the 
nonanalyticity of the coupling constant in all stages of the process for n > 1. 



IV. CONCLUDING REMARKS 

In this paper we have computed the self-energy of the neutral Xcf)^ model up to two- 
loop order at finite temperature. We have used MSCR in order to overcome the problems 
raised at finite temperature which demand resummation and renormalization. Thus, the 
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finite thermal mass has been consistently dressed up to two-loop order in the perturbative 
expansion and infinitely many other diagrams nonperturbatively. A comparison with other 



results in the recent literature may be done. In the second reference of the self-energy 
of a model having cubic and quartic interactions has been calculated at two-loop order. 
The model was used to illustrate the connection between multiloop self-energy diagrams 
and multiple scattering in a medium. The renormalization was discussed only for the first 
calculation of the self-energy i.e., with the vacuum mass m in the loops though the replace- 
ment of m by a thermal mass rriT in the loops is considered. In the 2PPI expansion 
has been used to compute the effective potential and the effective mass at two-loop in the 
model. Although they obtain renormalized quantities, there are mainly two differences. 
The first is that the corrections the 2PPI expansion take into account are not the complete 
two- loop corrections. The double scoop diagram is missing. As discussed in Ref. [0, in or- 
der to respect the symmetries of the Lagrangian, one must retain all diagrams to the given 
number of loops. The second difference is that they calculated the setting sun diagram at 
zero four-momentum. By definition this does not give the pole and necessarily the physical 
mass. Our work is more closely related with Parwani's resummed perturbative expansion 
P^ . He was concerned in obtaining consistently the pole of the effective propagator up 
to order g"^ (A^) taking into account temperature dependent counterterms as we did. The 
basic difference between his work and this one is that he used a mass corrected only by the 
one- loop diagram (which corresponds to our M|, Eq. (|35|)) in the calculation of the thermal 
mass up to two-loop. However, this assured the computation of the thermal mass up to 
order ^f^, as he intended. 

Although this model is very simple, the formalism applied here ought to be used in more 
realistic theories since the interaction and characteristics present here may find a place there. 
A next step could be the application of some properties of the MSCR in (more involved) 
gauge theories. 
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Figure Captions 

Figure 1: The one-loop contribution to the self-energy. 

Figure 2: The "extra" one- loop diagrams: (a) the mass counterterm diagram and (b) 
the vertex counterterm diagram. 

Figure 3: Two-loop contribution to the self-energy: (c) the "double scoop" diagram and 
(d) the "setting sun" diagram. 

Figure 4: A part of the complete self-energy diagrams in the superdaisy sum. The extra 
diagrams have been omitted. 



16 



extra 



n(M) = 






